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ABSTRACT. Tobacco smoking impairs CD4+ T-cells which coordinate the immune system. 
In this paper we develop a mathematical model to study how tobacco smoking affects the 
HIV infection of CD4+ T-cells. We derive the model equilibrium states and investigate 
conditions for their stability. Analysis shows that, disease free equilibrium is locally stable 
when Ror, < 1 which implies low smoking rate and when Ror, < 1 which imply low HIV 
infection rate, it is globally stable when Ror, < 1. Endemic equilibrium is asymptotically 
stable when there is high smoking rate in which Ror, > 1 and high HIV infection in which 
Ror, >1. 

The model is numerically simulated, results indicate that increasing smoking rate decreases 
uninfected T-cells and increases the number of HIV infected T-cells and free virus. The 
reduction of uninfected T-cells and the increase of HIV infected T-cells and free virus com- 
promises the immune system and expose a HIV smoker to opportunistic infections which 
marks a disease stage. To prolong and improve the quality of life for people living with HIV 


smoking cessation programs should be emphasized. 


1. INTRODUCTION 


The immune system is responsible for host defense against any external or internal ar- 
gent. In performing its defensive role, the immune system responds to a variety of antigens 
using T-cells with receptors for specific antigens and recognize body’s own antigens and 
activate harmful autoimmunity to the pathogens [10]. CD4+ T-cells in particular play a 
central role to protect the immune system through helping B cells to make antibodies, stir- 
ring macrophages’ ability to build up improved microbicidal activity, recruiting neutrophils, 
eosinophils and basophils to the location where there is infection or inflammation [38]. 


Infection of T-cells by HIV virus and impaired antigen response due to tobacco smoking 
have a remarkable negative effect to the immune system [35] and [37]. T-cells which are 
unresponsive due to tobacco smoke and are infected with HIV cannot perform their defense 
role because they cannot respond to the stimulation [19], [28] and have HIV DNA copy tran- 
scribed into their DNA [26]. Compounds such as nicotine, tar, carbon dioxide and carbon 
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monoxide from tobacco smoke induce anergy to T-cells and reduce their efficiency in re- 
sponding to body antigens and in producing cytokines and chemokines which are important 
for coordinating the immune system. Apart from reducing their defense potentials, T-cells 
are also used by HIV to further its reproduction during replication. T-cells which are af- 
fected by tobacco smoke and HIV fail to perform their normal functions and make the host 
vulnerable to various infections [27]. 


Although tobacco smoking has a devastating effect to the immune system, its effect on 
the HIV infection of T-cells is less known and has not been given attention whereas smok- 
ing rate of people who are living with HIV is increasing [2], [25], [36]. Mehta et al. [19], 
Sopori [27] and Kalra et al. [11] studied the effect of tobacco smoking on the immune 
system. Reviewing the effect of tobacco smoking, Mehta et al. [19] concentrated on the im- 
munosuppressive effects of tobacco smoke and mechanism through which tobacco smoking 
can compromise the immune system. On T-cells, which are the interest of this work, Mehta 
et al. [19] point out that smoking alters T-cells responses and reduce their defense potentials. 


Sopori [27] as well addressed the effects of cigarette smoke on the immune system. The 
study asserts that chronic exposure to tobacco smoke affects a wide range of immunological 
functions including responses from innate and adaptive immune systems. According to the 
study, nicotine and tar being the components of tobacco smoke are main cause of T-cells' 
insensitivity, these two compounds nicotine and tar also account for the immunosuppressive 
and addictive effect of tobacco smoking. Kalra et al. [11] investigated the effect of cigarette 
smoke on the immune responses. The study put forward that chronic exposure to cigarette 
smoke affects T-cells responses, this can account for diminished T-cells proliferation and 
antibody responses within the host. 


This study aims at using mathematical model to explore how tobacco smoking affects dy- 
namics of HIV in T-cells and how T-cells accommodate tobacco smoking and HIV. The 
recruitment rate of T-cells in our model is similar to the recruitment rate of T-cells in the 
model developed by Das et al. [5], however, tobacco smoking, productively and unproduc- 
tively T-cells differentiate the two models. Due to growing number of people who are living 
with HIV and smoke, this paper is interested to find out how tobacco smoking affects the 
HIV infection of CD4+ T-cells. 


2. MATERIALS AND METHODS 


2.1. Model development and Analysis. The model for the effect of tobacco smoking 
in the dynamics of HIV in T-cells divides T-cells into three classes, uninfected T-cells X, 
smoking impaired T-cells X; and HIV infected T-cells X2. The free virus population is 
represented by V. In this work, CD4+ T-cells which are the HIV target [35] are generally 
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referred to as T-cells. 





Uninfected T-cells are recruited at a rate i which is a decreasing function of the virus 


V [23] and [24], this is due to the ability of HIV virus to infect T-cells in the progenitor [5]. 
Uninfected T-cells acquire HIV and smoking impairment following a contact with free virus 
and smoking impaired T-cells at rates £1 and y respectively. 


Smoking impaired T-cells X, grow when uninfected T-cells acquire smoking impairment 
at a rate y, suffer HIV infection following a contact with free virus at a rate £5 and further 
suffer smoking induced mortality at a rate a. 


HIV infected T-cells X5 replenish at rates 6, and f» following a viral contact with unin- 
fected T-cells and smoking impaired T-cells respectively. The class of HIV infected T-cells 
X» also suffers HIV induced mortality at a rate jj. The parameter u represents T-cells 
natural mortality for all T-cells compartments. 


Free virus V increase when HIV infected T-cells replicate at a rate Nj; where N repre- 
sents number of free virus released by a single T-cell. Free virus diminish when they fuse 
with smoking impaired and uninfected T-cells at rates 91V X and G2V_X, respectively. How- 
ever, they suffer natural mortality at a rate ju. 


The model assumes that free virus prey on T-cells with 6ıV X and BV X, as functional 
responses respectively [5], tobacco smoking is treated as a disease and it is transmitted from 
one cell to another through contact; upon HIV infection, T-cells are assumed to be pro- 
ductively infected hence can produce free virus [26] and [29]; free virus are not affected by 
smoking; the contact between T-cells and free virus is through mass action principle. 

The state variables and parameters are defined in Tables 1, 2 respectively. The interaction 
between T-cells, smoking and free virus are shown by the Figure 1. 


TABLE 1. Variables description 
Variable Description 
X Uninfected CD44- T-cells 
Xi Smoking impaired CD4+ T-cells due to smoking 
X» HIV infected CD44- T-cells 
V Free virus 











J. IRUNDE, L. LUBOOBI, Y. NKANSAH-GYEKYE 


TABLE 2. Parameters’ description 





Parameter Description 





A constant for CD4+ T-cells recruitment 

k half saturation constant for the virus 

y rate of acquiring smoking impairment by CD4+ T-cells 
m natural mortality rate for CD4+ T-cells 

£1 rate of HIV infection for uninfected CD44- T-cells 

f» rate of HIV infection for impaired CD44- T-cells 

I HIV induced death rate for CD4+ T-cells 

[n natural mortality rate for free virus 

aq smoking induced death rate for impaired CD4+ T-cells 














pon 


FIGURE 1. Interaction of T-cells with smoking and free virus 
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Combining formulations and assumptions, the system that describes tobacco smoking 
effect on the HIV infection of T-cells is given by 


dX A 

(1a) dt = Ly = PVA — yXıX — uX, 
dX 

(1b) s = XX — BV X, — (a + u)X3, 
dX. 

(1c) Fre = BV X + BoV X — (u + ui) Xa, 
dV 

(1d) a = Nui X, — p1V X — 89V X, — lV. 


2.1.1. Invariant Region. To obtain the region where solutions of the model system (1) are 
feasible, we adopt the approach in [29] and the reference therein. The population for T-cells 
denoted by Xr and free virus V are considered where we have: 


(2) Xr= X +X +X. 


dXr 


ze < T — uXr, form = 


k 
The solution is given by: 
POE i («o — Ter 
H H 
for Xr(0) = X (0) + X1(0) + X2(0). 


Analysis of the solution X7 considers two cases; X7(0) > Z and Xr(0) < 
u 


for X7(0) > =; ee (xv) a z) ett and 
H Bu 
(3) 
forkle X. («o E z) ett < T, 
å n H H 


Since lim sup,;_,., («to — z) e"! — 0, 
(4) H 
limsup, sœ Xr < —. 
u 


From the definition of X7, it follows that: 
(5) lim sup X2 < E 
t—oo H 


Substitution of equation(5) in (1d) leads to: 


N 
(6) lim sup V < n 


t—>o Hy 
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Therefore we conclude that, the model system (1) is positive invariant in the region 
N 
(7) Dep V E ES Ce ay, re am) ay aa 
H mm 


Solutions for the model system (1) which are starting on the boundary of the region I 





enter the interior of the region and remain bounded. Results for existence, uniqueness and 
continuity for the system (1) hold in the invariant region I. Since the model is well posed we 
can therefore consider the flow generated by the model system (1). The result is summarized 
in the following theorem: 


Theorem 1. : All the solutions of the model system (1) enter the region 
N 
pert x Vier ake ee EP 
H "nm 


2.2. Disease free equilibrium and Reproduction number. We obtain disease free equi- 





librium when X, = X; = V = 0 and it is given by 
(8) Eo (X, X1, X3, V) E (5.0.0) i 
T 


In the absence of tobacco smoking and HIV, T-cells within human host assume their normal 
growth and function. The body at this state can fight against any pathogen that plans an 
invasion. At disease free equilibrium, the system of equations (1) reduces to 


— =n pX. 
(9) a 


The graph of T-cells showing asymptotic growth behaviour in the infection free state using 
parameters 7 = 20 and u = 0.02 [24] is illustrated in Figure 2. 

















i i i i i 
o 20 40 Ti 60 80 100 
Cursor position: (-5.69, ime 
The backward orbit from (-0.94, 3.824002) left the computation window. 








FIGURE 2. Behaviour of T-cells in the absence of smoking and free virus 
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2.2.1. Basic Reproduction number. To measure the effect of tobacco smoking in T-cells, the 
basic reproduction number which predicts the potentials of the disease [22] is computed . 
It determines how the two epidemics influence each other and which epidemic in particular 
accelerates the other. The basic reproduction number is defined as secondary infections 
caused by a single infected cell when introduced in the entirely susceptible population of 
cells [33]. To compute basic reproduction number Ror, the next generation method as 
described by van den Driessche and Watmough [31] is adopted. Using this approach, new 
infections and transfer terms are denoted by F; and V; respectively. The basic reproduction 
number Ror is given by: 





(10) Ror = p (FV"), 
that is the maximum eigenvalue of the matrix FV, where 
OF; OV; 
= —— (Xi = X 
F ax, | o) and V aX, | o) 


Using the model system (1) and the approach in van den Driessche and Watmough [31], F; 
and V; are defined as follows; 





yXı X 
(11) Fi = | AVX + 8:VXı 
Np X» 
and 
(at p1) Xi 
(12) Yi = | (u+ i) Xo 
HV 
whose eigenvalues A satisfy the equation 
(13) IFY! — AI| = 0, 
I is an identity matrix. The product of FV! is given 
or. 0 
wa + 1) 8 
T 
(14) FV = 0 0 f 
Hio 
Np 0 
H c pg 
Thus the basic reproduction number Ror is therefore given by 
(15) Ror — Max {Ror,, Ror, } " 
where 
YT Np Pir 
(16) Ror, = ————~ and Ron, = 4| ——————.. 
"plat p) 3 pps d ua) 
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The partial reproduction numbers Ror, and Rop, are representing new infections due to 
smoking and HIV respectively. Ror, depends on smoking impairment rate, the average time 
smoking impaired T-cell spends before death, life expectancy for T-cells and recruitment 
of T-cells in the absence of free virus, Rory, depends on the number of free virus which are 
produced by T-cells’ replication, HIV virus infection rate, life expectancy of T-cells and free 
virus and the recruitment of T-cells in the absence of free virus. 

We plot both partial reproduction numbers Ror, and Ror, against recruitment rate using 
parameters in Table 3 to compare the two epidemics in case of producing new infections. 
As indicated in Figure 3, tobacco smoking produces many new infections compared to HIV 
infection. 
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FIGURE 3. Comparison of smoking and HIV new infections 
2.3. Stability of Disease Free Equilibrium. 


2.3.1. Local Stability of Disease Free Equilibrium. The basic reproduction number Ror ob- 
tained by decomposition technique can be used to determine local stability of disease free 
equilibrium [8]. Disease free equilibrium is locally stable if Ror < 1 and locally unstable if 
Ror > 1. From the definition of Ror in (15), Ror < 1 will mean either Ror, < 1 or Ron < 1. 
Ror, < 1 if smoking rate is low and Ror, > 1 if smoking rate is high. In the same way 
Ror, < 1 if replication rate for T-cells is low and Ror, > 1 if replication rate is high. To 
analyze the local stability, we state and prove the following theorem: 


Theorem 2. : The disease free equilibrium for the model system (1) is locally asymptotically 
stable when Ror < 1 and is locally asymptotically unstable when Ror > 1 
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Proof 


To prove local stability of the model system (1), we linearize the system around a disease 
free equilibrium and obtain the following matrix: 


YT T Pır 
=i Ed 0 åt 
IE k p 
0 SS 0 0 
H 
(17) Jeo Z 
0 0 —( + u) Ee 
fn 
0 0 Nie . dn. 
H 


For disease free equilibrium to be locally asymptotically stable, eigenvalues of Jacobian 
matrix should be negative. From the first and second columns of the Jacobian matrix we 
obtain —y and D we u which we re-write as; 

u 





(18) e+) (7 1) = e (Ron - 1). 


It is negative if and only if 


(19) Ror, « 1. 
The Jacobian matrix (17) is now a 2 x 2 matrix; 
T 
" -m+ BR 
20 Jg, = 
1 T 
Np -Ar — py 


To establish local stability of disease free equilibrium, trace and determinant are computed. 
Local stability holds if the trace is negative and determinant is positive. Trace of the matrix 
Jp, is given by 














(21) Jm =~ p) -ÉT py <0. 
The determinant of the matrix Jp, is: 
N 
(22) det Dp, = PTEE 4 pp 4 pm) - SHA, 
H H 
which simplifies to 
dn iU Hon. 
(23) Du 
= T (1 + Ror,)(1 NE Ror). 


Holt 
It follows that 


(24) det Jg, > 0, iff Ror, «1 
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From equations (19) and (24), disease free equilibrium is locally asymptotically stable when 
Ror, < 1 and Ror, < 1. This completes the proof. 


2.3.2. Global stability for Disease Free Equilibrium. Local stability of disease free equilibrium 
gives the overview for the behaviour of the disease around the disease free equilibrium point. 
To provide comprehensive understanding for the behaviour of the disease beyond disease free 
equilibrium point, global stability analysis for disease free equilibrium point is carried out 
and we adopt the method in Castillo-Chavez et al. [3]. Using this method, the system (1) is 


re-written as 








dXn 
= C (Xn — Xae) + C1X;, 
e d, CX 
dt > 24°45 


where X,, represents non-transmitting classes, X; represents transmitting classes and X. a 
disease free equilibrium point respectively. C, Cı and C are matrices to be obtained from 
the system (25). Global stability is guaranteed if eigenvalues of matrix C are negative and 
C» is a Metzler matrix [20] defined mathematically as C»(r;j) > 0 Vi # j. Adopting the 
form in equation (25), the system (1) is written as; 


A 7 = 
(26) E MAP Rot de -c(X-z)«a X; 
k+V H 
V 
and 
yXı X = 83V X4 = (a + py) Xi X1 
(27) BV X + 83V X4 — (ji + pp) Xy =C | Xo 
N mX — BV X — B3V X1 — jV V 
Matrix C is 1 x 1 matrix and it is C = —y with eigenvalue —44. Matrix Cj is a 1 x 3 matrix 
given by 
C= (= RE M ar) 
h k p 
Matrix C5 is also given by 
a — (a+ u) 0 0 
T 
(28) Cy = 0 —(u + Hu) " 
pn 
0 d -T - Hv 
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where d = Np. The matrix (28) has non-negative off diagonal elements and all diagonal 


T 
elements are negative except EE (a + u) which is further written as 
u 


Tf (adele a) (Ren e d ys 
u 


(a + u) (Ror, — 1) < 0, iff Ror, < 1. 


(29) 


Ror, < 1 implies low smoking rate and Ror, > 1 high smoking rate. From this findings 
it is seen that low smoking rate guarantees global stability of disease free equilibrium and 
high smoking rate renders global instability of disease free equilibrium. These findings are 


summarized in the following theorem: 


Theorem 3. : The disease free equilibrium is globally asymptotically stable if Ror, < 1 and 
it is globally unstable if Ror, > 1. 


2.4. Endemic Equilibrium. To determine how the two epidemics interact within the host, 
we compute the endemic equilibrium and derive the condition for its stability. To compute 
endemic equilibrium, system (1) is written in the following form: 


30a) eO ea cades): 


l krV 

(30b) YX X — B3V X4 — (a + u)Xı = 0, 
(30c) BV X + 83V X1 — (wt po) X5 = 0, 
(30d) Nytt Xo — BV X — BV X1 — jV — 0. 


Beginning with equation (30b), we have 


yXiX — bV X; — (at u) X1 = Xi [1 X — bV — (a + u)| = 0, 
(31) B3V* + uy 


AG =0, X* SS aa Soe. 
Sr 
Xj = 0 represents disease free equilibrium when Vj* = 0. To obtain XT we substitute X* in 
(30a) where the equation now becomes 


A BoV* + Uy BoV* + Uy 

ss TANN sj PE ENES 

ae Bi ( 7 ) 1 (b2 V* + u) Å F ) , 
NA A | Pov" tu 


Ai (k + V*)(82V* + u) yo? 


In a simplified form, we have 





Ay — (k + V*)(BoV* + u1)(B1V* + n) 
qlk + V*)(BoV* + u) l 

From equation (30c), we can obtain Xž as follows; 

| BiV*X* + BoV*XT 


U2 


(32) X: = 





(33) x 





uo = U + H. 
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Substitution of X* and X7 into equation (33) results into the expression for X5 
(k + V*)(BV* + u1) + (181 — n82)V* + Ay BV" 
uzy(k + V*)(83V* + ui) 
All variables, X*, Xf and X3 are in terms of V*. To find V*, X*, Xt and X5 are substituted 
in equation (30d) where we obtain 
d(k + V*)(B3V* + ui)(ui Bi — n82)V* + AyBoV* — BYV"(BoV* + ux) 


(35) uoy(k + V*)(85V* + u1) y 
, BV" (k + V*)(AYV" + p) (b2 V* + u1) — Agp V* diee 


ylk + V*)(BoV* + u1) 


Expansion and rearrangements of terms in equation (35) yields a third degree polynomial in 





(34) Xi = 








variable V*; 
(36) aV* +bV** + cV* — 0, 


written as 

V*[aV*? + bV* +c] = 0, 

a = By[(d — u3)(u181 — 182) — Yuva], 

b = (kB + ui)|(ui£1 — u82)(d — u2) — usua], 

c = (d — u3)|y AB» + uik(uiBy — uf2)] — ypsuiusk. 


V* = 0 represents disease free equilibrium when Xf = 0. The equation 


(37) 


(38) aV + bV* - c — 0, 


represents endemic equilibrium. Using general quadratic formula in equation (38), 


Ey vb? — 4ac 


2a 
— 2 
V* > gig ete å 


2a 
—b + Vb? — 4ac 
a 


5 > 0 iff ac < 0. 
ac<O0ifa<Oorc<0. 


V 











YHyU2 
(d = uz) (u1 bı = 2) 


a < 0 iff > lor 





(39) 
HyUu2 = AB» 
(uiBy — up2)(d — u3) ^ uyk(u8y — ufa) 


From equation (39), a < 0 if there is high smoking impairment rate y and low HIV infection 





c « D iff y ou 


rate £5 for smoking impaired T-cells and c < 0 if there is high smoking impairment rate y 
and high HIV infection rate £» for smoking impaired T-cells. Therefore either a < 0 orc < 0 
but not both. Endemic equilibrium is defined by the set I'* such that 


(40) pex. XV 
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where the values of X*, X1 and X5 in terms of V* are given in equations (31), (32) and (34). 
This result is summarized in the following theorem: 


Theorem 4. : High smoking impairment rate provides the necessary condition for endemic 
equilibrium I* = {X*, Xt, X3,V*} of model system (1) to exist. 


2.4.1. Global stability of endemic equilibrium. Stability analysis explores the behaviour of 
the epidemic near the equilibrium points. Solutions which start near the equilibrium point 
and remain near for all times are stable solutions and represent a stable behaviour of the 
epidemic. Solutions which start near the equilibrium point and converge to the equilibrium 
point are asymptotically stable and they represent asymptotically stable behaviour of the 
epidemic. The solution which show neither of the two behaviours represent unstable be- 
haviour of the epidemic. 


The Lyapunov method and LaSalle’s invariance principle have been widely used to investi- 
gate the stability of epidemics. Using Lyapunov method and LaSalle’s invariance principle 
different functions are constructed depending on the types of the models formulated. 


Logarithmic Lyapunov function [7] was developed to analyze stability of Lotka-Volterra 
systems. Korobeinikov [15] used it to analyze stability for SIS, SIR and SIRS models. 
Vargas-De-Leon [32] used composite quadratic Lyapunov function to prove global stability 
for SIS, SIR and SIRS models and Korobeinikov [13] and [14] used explicit Lyapunov 
function to analyze SEIR and SEIS epidemic models. 

To analyze the endemic equilibrium for the model system (1), we adopt explicit Lyapunov 


function 
(41) H(yi) = 3 wilyi — yi In yi) 
i=1 


where w; are constants that are to be carefully selected, 4; represents number of cells in class 
i and y; represents equilibrium state in class i. According to the approach in McCluskey 
[18], using the model system (1) we now define the explicit Lyapunov function 


Ht XS XV) ET KMV SOR by: 


A(X, X41, Xo, V) = w(X — X* In X) + w(Xı — Xi In X1) 


(42) 
+w3(Xo = X35 In X5) + w4(V — V*ln V). 


We assume that for i = 1,...,4 w; > 0, function H and all w; are continuous and differentiable 
in I. The value of the function H in I" is zero. 
Differentiating the explicit Lyapunov function H with respect to time leads to the following 
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equation 
dH | X*\d Xi Y dX, X35 X dX» 
(43) dt m (1-3) Fre 1-2) G uy (1-32) Te 
2 
Vj dt’ 


which further gives; 
dH 
dt 


a 


+ [XX = bV X4 = (a + iu) X4| 





(44) 
+ [63V X + 89V X4 — (u + por) X3] 








Tr p E [Np X» — BV X — BV X1 — p, V] , 


At endemic equilibrium, equation (44) becomes; 




















dH X* 
FE (i- x) [01V*X* + Xt X* + uX* — BV X —yX4X — uX] 
AX 
+w | 1 — X, [b2V*X{ — (a + u) XT — BV X1 — (a + u)Xıl] 
45 X3 x 
(45) +w; | 1— X, (u 1) X3 — (u pr) X2] 
V* 
+w4 | I — v) [61V * X* + B9V* Xt + u,V* — 8V* X* — B9V* Xi 
=p V*]. 
Rearrangement of terms and further simplification gives 
dH n | Qa - X? 
WM LES LIED ERE 2 Fu) X : 
Xə — X> V- V* 
—wa(p 4 i) oe 2) mm ee 
X — XY (VX —V*x* X — X*) XX — XIX" 
(46) “ang XX ) mn XX) GX - XIX") 








(x - X (bx, VX) (v -vtvx - vx") 

















—wW2/B2 X, are WAP1 V 
(V —V*)(VX, —V*XI1) 
—w4 p2 , 
V 
which is written as 
dH DOE bis b OR dn 
Es EE 
2/312! 2 = 
— j — T F(T 
w(u + pi) X; 4H V +F(T) 


where 


r= {(X, X, X2, V) > 0j 
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and 


je eu LI E UA 


X 
(X; — X1) (VX, - V* X1) 








(V — V*) (V — V*X*") 
V 





— wap 





—w b2 





X 
(V —V*)(VX; —V*XI) 
V i 
According to McCluskey [18] and Korobeinikov [14] F(T) is non-positive in I and therefore 


dH dH 
F(T) € 0 for all T. The derivative — < 0 in I, equality holds when I = I*. Since — < 0 


—w4p2 


dt dt 
H dH 
for all I and eH = 0 when I = I, it means that the largest invariant set in IT when pi 0 


is a singleton (I*) which is the endemic equilibrium. Hence by LaSalle’s invariance principle 
[16] and [17] it means that the endemic equilibrium point I* is asymptotically stable in I 
when Ror > 1. This result is summarized in the Theorem 5: 


Theorem 5. : Endemic equilibrium of the model system (1) is globally asymptotically stable 
when Ror > 1 and it is globally asymptotically unstable otherwise. 


3. NUMERICAL SIMULATIONS 


To accomplish the findings of this study, numerical simulations are carried out to analyze 
the long term behaviour of the two epidemics. Analytical results show that tobacco smoking 
has devastating effect to HIV infection of T-cells. High smoking rate provides a condition for 
endemic equilibrium to exist. Stability analysis of disease free equilibrium shows that tobacco 
smoking provides necessary condition for local and global stability of disease free equilibrium. 


To study further the effect of tobacco smoking on the HIV infection of T-cells, the gen- 
eral dynamics is examined to reveal the features for the epidemics in T-cells. Then the 
specific features which give the overall behaviour for the two epidemics are highlighted. 


The life expectancy of a HIV infected patient according to Nakagawa et al.[21], is 9-10 
years from primary infection. We will use this time span in our simulations. However, to 
critically study provide a clear picture of the two epidemics, we will increase the time span 
to 20-30 years. The parameter values from the literature are used. According to Kirshner 
and Webb (1997), the normal recruitment of CD4+ T-cells is 10/mm? per day as also seen 
in [29] and [34]. However, due to HIV infection this rate falls from 10/mm? per day to 
3/mm? per day [12]. In this work the recruitment of 3/mm? per day which is for a person 
living with HIV [12] is used. Basing on the bahaviour of tobacco smoking and HIV epidemic, 
estimated parameter values will also be used due to limited studies which have addressed 
tobacco smoking and HIV. The parameters are listed in the Table 3. 
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TABLE 3. Parameter Values 








Parameter Description Value Source 
A constant for CD4+ T-cells recruitment 1095 year” 12] 
k Half saturation constant 12 5 
Bi HIV infection rate for T-cells 0.000292 — p 
virus year 
y Smoking impairment rate for T-cells 0.0025 year-! 9 
u natural mortality rate for T-cells 0.02 year! 35] 
Bo HIV iotesion ate Gp Pade. COS 
virus year 

I HIV induced mortality 0.0245 year! 9 
iin natural mortality rate for virus 10 year! 24] 
q Smoking induced mortality 0.00024 year! 9 











The general dynamics of HIV in T-cells in the presence of tobacco smoking is plotted in 
Figure 4. It reveals that smoking impairment grows faster within the first five years. HIV 
infection raises when smoking impairment is on the decline. At this point, the uninfected 
T-cells decline dramatically and reach their lowest value as portrayed by Figure 4. 
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FIGURE 4. Variation of T-cells in the presence of smoking and HIV 

The growth of HIV infection which results into increase of HIV infected T-cells while 
smoking impaired T-cells are declining shows that smoking impaired T-cells are infected by 
HIV at a high rate. The graph of smoking impaired T-cells against HIV infected T-cells 
shows a corresponding increase of HIV infected T-cells with respect to smoking impaired 
T-cells. Initially smoking impaired T-cells increase but it reaches a point where it ceases. 
At this point HIV infected T-cells grow to their maximum and decline thereafter to attain 
the value that does not change as illustrated in Figure 5. However, when time span tripled 
HIV infected T-cells decline to their lowest value due to replication and natural mortality. 
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FIGURE 5. Variation of Smoking impaired and HIV infected T-cells with re- 
spect time and Variation of smoking impaired T-cells with HIV infected T-cells 
Smoking impairment rate is then varied to see its effect on the dynamics of HIV in T-cells. 

Our aim in this case is to observe how does this variation affects uninfected T-cells, smoking 
impaired T-cells and HIV infected T-cells. We also look at the relationship between impaired 
T-cells and HIV infected T-cells and the corresponding effect on free virus. 
The uninfected T-cells decrease as smoking impairment rate increases, free virus also increase 
per smoking impairment increase as demonstrated in Figure 6. Decreasing of the uninfected 
T-cells and increasing of free virus can worsen T-cells’ profile and pave a way for various 
infections which can accelerate HIV infection to the disease stage AIDS and ultimately death 
of the human host. 
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FIGURE 6. Variation of uninfected T-cells and free virus with smoking im- 
pairment rate. 

'The smoking impaired compartment is also studied when smoking impairment rate varies. 
This variation results into a corresponding increase of smoking impaired T-cells, this is 
observed within 10 years. Tripling time span results into increase of impaired T-cells but as 
smoking impaired T-cells increase they take less time to decline to their minimum. This is 
because the more the T-cells are impaired by smoking the faster they acquire HIV infection 
and this is the reason why they diminish for a short time as depicted in Figure 7. The 
result of high HIV infection rate to smoking impaired T-cells is reflected in HIV infected 
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T-cells because as smoking impairment rate increases there is a corresponding increase in 
HIV infected T-cells as demonstrated by Figures 8 and 9. 
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FIGURE 7. Variation of impaired T-cells with respect to smoking impairment rate. 



































































800 800 
y=0.0002 
y-0.004 
700r 700 y=0.0008 
600r 600 
å 500r 2 500 
E å 
E 400 - B 400 
= 300} = 300 
rI 
200r 200 
100r 100 
0 o 
o 5 10 15 
Time[years] Time[years] 
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FIGURE 9. Variation of HIV infected T-cells with respect to smoking impaired 
T-cells. 


4. DISCUSSION AND CONCLUDING REMARKS 


A simple model for tobacco smoking effect on the HIV infection of T-cells is formulated, 


equilibrium states and their condition for stability are derived and basic reproduction number 
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is computed. At disease free equilibrium, T-cells grow asymptotically. The basic reproduc- 
tion number which is given as the maximum of partial reproductive numbers Ror, and Ror, 
shows that among the two epidemics tobacco smoking produces many new infections in T- 
cells compared to HIV infection as demonstrated by the Figure 3. 


Analytical result shows that, equilibrium states (disease free and endemic equilibriums) 
exist and smoking affects the existence of endemic equilibrium. Stability analysis of the 
equilibrium states shows that; when Rory, < 1 and Ror, < 1 disease free equilibrium is lo- 
cally asymptotically stable and it is globally stable when there is low smoking rate for which 
Ror, < 1. When Ror, > 1 which implies there is high smoking impairment rate and Ron, > 1 
which implies there is high HIV infection rate, the endemic equilibrium is globally stable. 


The model is numerically simulated to study the effect of tobacco smoking. We observe 
that as tobacco smoking impairment rate increases, uninfected T-cells reduces dramatically, 
HIV infected T-cells and free virus increase significantly. Decreasing of uninfected T-cells 
is due to smoking impairment and HIV infection. HIV infection of uninfected T-cells and 
smoking impaired T-cells results into many HIV infected T-cells which replicate to release 
free virus. As many T-cells are impaired by tobacco smoking, HIV infected T-cells increase 
correspondingly. This is due to the fact that smoking induced anergy affects the normal 
functioning of T-cells by making them unresponsive and easily attacked by HIV. This is 
reflected by increase of HIV infected T-cells when smoking impairment rate increases as 
depicted by the Figure 8. 


Reduction of uninfected T-cells and increase of HIV infected T-cells and free virus present a 
compromised immune system for which a host is exposed to opportunistic infections [6] such 
as tuberculosis, oral candidiasis [4] and cancer. Tobacco smoking can therefore promote HIV 
infections in T-cells and deteriorate the T-cells’ profile as proposed by Valiathan et al.[30]. 
For individuals living with HIV to live longer and improve the quality of life, the study calls 
for smoking cessation programs to emphasize the importance of abstinence from tobacco 
smoking and other risk behaviours. 
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